The slow-mode nature of compressible wave power in solar wind turbulence 
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We use a large, statistical set of measurements from the Wind spacecraft at 1 AU, and supporting 
synthetic spacecraft data based on kinetic plasma theory, to show that the compressible component 
of inertial range solar wind turbulence is primarily in the kinetic slow mode. The zero-lag cross 
correlation C(Sn, SBu) between proton density fluctuations Sn and the field-aligned (compressible) 
component of the magnetic field <5B|| is negative and close to —1. The typical dependence of 
C(5n, 8B\\ ) on the ion plasma beta pi is consistent with a spectrum of compressible wave energy 
that is almost entirely in the kinetic slow mode. This has important implications for both the nature 
of the density fluctuation spectrum and for the cascade of kinetic turbulence to short wavelengths, 
favoring evolution to the kinetic Alfven wave mode rather than the (fast) whistler mode. 
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Introduction. — The inertial range of solar wind tur- 
bulence is comprised of a mixture of incompressible and 
compressible motions, with at least 90% of the energy due 
to the incompressible component [l[ . If these fluctuations 
are interpreted as some mixture of the three MHD linear 
wave modes, then Alfven waves are the dominant incom- 
pressible component, while slow and fast MHD waves 
make up the compressible component. These modes are 
distinguished by the correlation between the density and 
parallel magnetic field fluctuations: fast waves are pos- 
itively correlated, slow waves are negatively correlated, 
and the density and parallel magnetic field fluctuations 
are both zero for Alfven waves. As the wave amplitude is 
increased to nonlinear levels, even to the limit that they 
form discontinuities or shocks, these qualitative proper- 
ties persist, corresponding to tangential and rotational 
discontinuities or fast and slow shocks Q. The MHD 
limit of strong collisionality, however, is not valid in the 
solar wind; therefore, collisionless kinetic theory is neces- 
sary to determine the properties of the wave modes. Each 
of the kinetic versions of the MHD linear wave modes, 
determined using the Vlasov-Maxwell linear dispersion 
relation, retain the qualitative correlations between the 
density and parallel magnetic field fluctuation described 
above jjjjj. In addition, these kinetic counterparts to the 
compressible modes may suffer damping from collision- 
less mechanisms 

Compressible fluctuations at inertial range scales A in 
the solar wind (1CP 4 Hz < f sc ~ v sw /X < 1 Hz; v sw 
is the solar wind speed, f sc is the Doppler-shifted fre- 
quency in the spacecraft frame) have been studied exten- 
sively, often interpreted as a mix of magnetoacoustic (fast 
MHD) waves and pressure-balanced structures (PBSs) 
0, H| ■ A PBS was first observed as an anti-correlation 
of thermal pressure and magnetic pressure at timescales 
of 1 h Q, and subsequent investigations found a sim- 
ilar anti-correlation between the density and magnetic 



field magnitude 0,@|. Theoretical studies of compressible 
MHD fluctuations in the low-Mach number, high-/? limit 
interpreted these anti-correlated density-magnetic field 
strength observations as nonpropagating "pseudosound" 
density fluctuations A more comprehensive investi- 
gation confirmed the general density-magnetic field anti- 
correlation, but also identified a few positively correlated 
intervals consistent with the magnetosonic (fast MHD) 
wave [Icj . Analysis of Ulysses observations found evi- 
dence for PBSs at inertial range scales in the high latitude 
solar wind (llTfl3j . Studies of the electron density up to 
f sc = 2.5 Hz also found pressure balanced structures, 
but interpreted these as ion acoustic (slow MHD) waves, 
and recognized that PBSs are simply the ion acoustic 
(slow MHD) wave in the perpendicular wavevector limit 
14j . Recently, measurements of the anti-correlation be- 
tween electron density and magnetic field strength indi- 
cated the existence of PBSs over timescales ranging from 
10 3 s down to 10 s 0. 

This Letter demonstrates that the compressible fluc- 
tuations in the inertial range are almost entirely kinetic 
slow wave fluctuations, suggesting that little turbulent 
energy is transferred from large scales to whistler fluctu- 
ations below the ion gyroscale. First, we show that the 
density-magnetic field cross-correlation C(Sn, 5B\\ ) in the 
solar wind is ~ — 1 and increases slightly with ion plasma 
beta, j3i . Then we demonstrate excellent agreement with 
synthetic (eigenfunction) data in which less than 10% of 
the compressible energy is due to fast waves. 

Measurements. — We use measurements from the Mag- 
netic Field Investigation (MFI) [l6[ and the Three Di- 
mensional Plasma (3DP) experiment [13] on the Wind 
spacecraft, in the unperturbed solar wind at 1 AU, dur- 
ing the years 1994-2004. The magnetic field is sampled 
at either 11 or 22 vectors/s (depending on the space- 
craft telemetry rate) then averaged down to the space- 
craft spin period (3 s). The ion moments are computed 
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FIG. 1. Example of waveforms of proton density fluctu- 
ations Sn (black with red dots) and parallel magnetic field 
<5B|I (green) fluctuations. This interval has (<5n) r ~ 0.6 cm -3 , 
0.4 nT, and C(Sn,5B\\) ~ -0.97, however, even in- 
tervals with smaller density fluctuations exhibit significant 
ant i-correlations. 



on board the spacecraft at 3 s cadence; protons are sep- 
arated from alpha particles by a fixed energy interval 
(which is occasionally adjusted in flight software) during 
the moment calculation; the solar wind alpha particle 
abundance is typically 3-5% of the proton number den- 
sity. We select 1,089,491 300-s intervals of ambient solar 
wind data (corresponding to spatial intervals of approx- 
imately L ps 450 km/s x 300 s = 135,000 km ~ 1350 
Pi, where pt is the ion Larmor radius) and the data is 
decimated by a factor of 10 (to 30 s cadence). There- 
fore, our data correspond to inertial range scales of ap- 
proximately kpi £ (5 x 10~ 3 ,5 x 10~ 2 ). After selection, 
the magnetic field data are averaged in 100-s windows to 
compute the local mean field Bo- The fluctuation field 
(5B is created by subtracting Bo and then rotated to a 
field-aligned coordinate (FAC) system defined by the Bo 
direction. In this new system, there is a compressible 
field fluctuation 8B\\ and shear components 8B± t i and 
5B± t 2- The most probable amplitude of the shear com- 
ponent SB± = {8B\ 1 + 8B\ 2 ) 1 ^ 2 is approximately 3.4 
times greater than most probable 8B\\ and corresponds 
to the Alfvcnic component of the turbulence. 

Proton density data 8n = n — uq is detrended over the 
same time intervals. Proton density is an integral over 
the (3 s) distribution function /(v), computed from par- 
ticle flux measurements in 16 individual energy channels 
[l7j . Since the counts are digitized discretely and the 
energy channels are also discrete (with AE/E ~ 0.2), 
the 3DP proton density moments have a finite dynamic 
range. To assess this, we evaluated the joint probabil- 
ity distribution (not shown) of (Sn) r and (8Bu) r , the 
RMS values. Below values of (5n) r ps 0.5 cm -3 , the joint 
pdf reverts from being well-correlated to a broader set 
of values. Furthermore, the histogram of (Sn) r alone 
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FIG. 2. Histogram of the cross-correlations C(8n, 8Bu) (grey, 
1,089,491 points total) and those above the measurement 
noise threshold Sn > 0.5 cm (blue, 119,512 points total). 



shows clearly an artificial (non-Poisson) lower cutoff at 
around this value, a cutoff not seen in (8B\\) r . We take 
(<Jn) r PS 0.5 cm -3 as the "noise" level of the density fluctu- 
ation measurement and restrict our analysis to intervals 
in which (8n) r > 0.5 cm -3 . This restriction reduces the 
dataset from 1,089,491 to 119,512 data intervals and bi- 
ases our sample to higher absolute densities. The most 
probable values of absolute density are ~4 cm -3 for the 
full and ~12 cm -3 for the thresholded dataset. The dis- 
tribution of plasma Pi is unaffected by the thresholding. 

We compute the normalized, zero-lag cross-correlation 
C(8n,8B\\) = (SnSBu)/ (8n)r(8B«) r , which has a range 
from —1 to 1. As described qualitatively above, we ex- 
pect that C(8n,8B\\) will be negative (positive) for slow- 
(fast-) mode MHD fluctuations (see Figure 1). Figure 2 
shows the distribution of C(8n,8B\\), both for all of the 
data and for the restricted dataset that exceeds the den- 
sity noise threshold (8n) r > 0.5 cm -3 . While the en- 
tire dataset peaks below (at C(8n,8Bn) ~ -0.5), the 
data with well resolved density amplitude levels peaks at 
C(8n, SBn) ~ -0.9. Figure 3 shows the joint histogram of 
C(8n,8B\\) vs. ion plasma beta /3,. The top panel shows 
the distribution of points, with a histogram of pi over- 
plotted. The middle panel is the joint histogram normal- 
ized to number of Pi points in each Pi bin. This shows 
clearly that C{8n,8B\\) is near —1 over the entire inter- 
val, increasing slightly with Pi to ~ —0.7 at pi = 10. In 
the bottom panel, the cumulative distribution shows that 
fewer than 10% of the intervals have C(8n, SBh) > 0. 

Synthetic Data. — A cubic synthetic plasma volume 
spanning scales 3 x 10~ 3 < kpi < 4.8 xl0~ 2 is constructed 
using a 32 3 grid. A spectrum of linear waves, with 90% 
of the energy in Alfven waves and the remaining 10% in 
a mixture of kinetic fast and slow waves, consistent with 
the observed fc _5//3 one-dimensional energy spectrum of 
the magnetic field fluctuations |£B|, is created in the vol- 
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FIG. 3. (Upper) Total distribution of the C(Sn,SB^) cross- 
correlation as a function of ion plasma beta /3i. The count 
in each /3; bin is overplotted (with the scale on the right). 
(Middle) Joint distribution of the C(5n, 5B^ ) cross-correlation 
normalized within each f3i bin. Black dots are the peak values 
in each j3i bin. (Lower) Cumulative distribution of C(Sn, SB^) 
with contours at 90%, 50%, and 10%. These distributions 
consist of 119,512 independent data intervals. 



ume using the linear eigenfunctions for these modes from 
the Vlasov-Maxwell linear dispersion relation 18 , 19[ ; see 
Klein et al. 2011 Q for more details. A fully-ionized 
proton and electron plasma is assumed, with isotropic 
Maxwellian velocity distributions, a realistic mass ra- 
tio m,/m e = 1836, equal ion and electron temperatures 
Ti = T e , and non-relativistic conditions vu/c = 10~ 4 . 
Taking the MHD limit kpi <C 1, under these conditions 
the normalized linear Vlasov-Maxwell eigenfrcquency de- 
pends on only two parameters, ui/{kvA) = uJ(/3i,#), the 
ion plasma beta f3i and the angle 6 between the wavevec- 
tor and the mean magnetic field 0]. Once (3i has been 
chosen, one needs only to specify the distribution of en- 
ergy in wavevector space. Compressible MHD turbu- 
lence simulations generate an isotropic distribution of 
fast waves and critically balanced distributions of Alfven 
and slow waves [2oT ]. Therefore, we initialize the fast wave 
energy isotropically, while the Alfven and slow wave en- 
ergy mimics a critically balanced distribution by setting 



all modes with /c|| > k^k 2 ^ 3 to zero, where fco corre- 
sponds to the scale of the plasma volume. 

Time series of density and parallel magnetic field fluc- 
tuations are constructed by sampling the synthetic data 
at a probe moving through the volume at an oblique angle 
with respect to the mean field (tests have confirmed in- 
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FIG. 4. Comparison of measured values of the C(5n,5B«) 
cross-correlation (black dots with FWHM error bars) to the 
synthetic data predictions for the ratio of kinetic fast wave 
to total compressible energy F. Best agreement is with 
F = 0.00, indicating that the compressible component of solar 
wind turbulence is almost entirely in the kinetic slow mode. 



sensitivity to the choice of angle). We then compute the 
cross-correlation C(Sn, 6Bn ) as above. Figure 4 shows 
C(Sn, SBh ) for several values of the ratio of fast wave en- 
ergy to total compressible energy F vs. ion plasma /3j. 
Peak histogram values (and FWHM error bars) are also 
shown. The solar wind data are in striking agreement 
with the synthetic data F = 0.00 curve. 

Discussion. — Figure 4 shows that the observed corre- 
lation is consistent with a statistically negligible kinetic 
fast wave energy contribution for the large sample used in 
this study. Note, however, that a very small fraction of 
the data intervals have a positive cross-correlation (see 
Fig. [5]), possibly indicating the presence of kinetic fast 
waves in these intervals. 

Previous analyses have generally dismissed the pos- 
sibility of kinetic slow waves because, in an isotropic 
Maxwellian plasma with warm ions, the collisionless 
damping via free-streaming along the magnetic field is 
strong [4|. However, in the limit k± 3> fc|| applicable to a 
critically balanced power distribution, the damping rate 
of the slow waves is proportional to the parallel compo- 
nent of the wavevector, 7 oc fc|| [l9j| . For exactly perpen- 
dicular wavevectors, the damping rate drops to zero — 
this perpendicular limit of the slow wave corresponds to 
an undamped, non-propagating pressure-balanced struc- 
ture (PBS). In compressible, strong MHD turbulence, it 
has been shown that the slow modes are cascaded pas- 
sively by the Alfvenic turbulence [2l|, [22J , so the energy 
cascade rate is related not to the slow wave frequency, 
but to the Alfven wave frequency. Therefore, the more 
nearly perpendicular slow waves (possibly with k^ well 



vg^/c^ 3 ) may be cascaded 
to smaller scales on the timescale of the Alfvenic fur- 



below critical balance, fc|| <C fc. 
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FIG. 5. Plot of frequency aj/fi; vs. wavenumber kpi for the 
collisionless versions of the fast MHD (red), Alfven (blue), 
and slow MHD (green) waves determined using the Vlasov- 
Maxwell linear dispersion relation. The range of scales con- 
sidered in this study is indicated. Our measurements suggest 
that the solar wind spectrum lies on the blue-green curves. 



bulence, while the collisionless damping of these modes 
remains weak. 

We offer the following physical model of the compress- 
ible fluctuations in solar wind turbulence. At inertial 
range scales kpi < 0.1, the density and parallel magnetic 
field fluctuations arise mainly from the kinetic counter- 
parts of the fast and/or slow MHD waves. The measured 
C(Sn, (5B|| ) cross-correlation at these scales demonstrates 
that the compressible fluctuations are statistically dom- 
inated by kinetic slow mode fluctuations, and the distri- 
bution of power in wavevector space of these slow modes 
mimics the critically balanced distribution expected of 
Alfvenic fluctuations 0] . These kinetic slow wave fluctu- 
ations may be cascaded as passive fluctuations to smaller 

| , and so are pre- 



scales by the Alfvenic turbulence [21 



dieted to exist down to the scale of the ion Larmor radius, 
and perhaps to even smaller scales. Thus, the evidence 
for PBSs over a range of timescales from 10 3 s to 10 s 



15j is explained by the presence of a distribution of ki- 
netic slow waves that is undergoing a turbulent cascade to 
smaller scales driven by the Alfvenic turbulence. In addi- 
tion, these kinetic slow modes suffer collisionless damping 
at a rate y sx kn , meaning that the more perpendicular 
the wavevector, the slower the damping rate. 

The lack of a statistically significant fast wave compo- 
nent in the inertial range of solar wind turbulence has 
implications for the cascade of energy to small scales. 
Because the nonlinear energy transfer is believed to be 
dominated by local interactions in wavenumber space, 
significant nonlinear energy transfer occurs only between 
waves with similar linear frequencies. Fig. [5] shows the 
real linear frequencies lu of the collisionless counterparts 
of the MHD fast, Alfven, and slow waves as a function 
of kpi. For the isotropically distributed fast waves (red), 
we plot the parallel (solid), 45° (dotted), and perpen- 



dicular (dashed) increase of the wavevector; slow (green) 
and Alfven (blue) waves follow the critically balanced 
path fc|| = k^k 2 ^ 3 , with the isotropic driving scale 
k pi = 10~ 4 for all cases. Since only the fast wave tur- 
bulent cascade is expected to nonlinearly transfer energy 
to whistler waves at kpi > 1, our analysis suggests that 
there is little or no transfer of large scale turbulent en- 
ergy through the inertial range down to whistler waves 
at small scales. 
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